Phase diagram of the D3/D5 system in a magnetic field and a BKT transition 
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We study the full temperature and chemical potential dependence of the D3/D5 2+1 dimensional 
theory in the presence of a magnetic field. The theory displays separate transitions associated with 
chiral symmetry breaking and melting of the bound states. We display the phase diagram which 
has areas with first and second order transitions meeting at two critical points similar to that of 
the D3/D7 system. In addition there is the recently reported BKT transition at zero temperature 
leading to distinct structure at low temperatures. 
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I. INTRODUCTION 

There has been recent interest in holographic descrip- 
tions of the phase structure of gauge theories in the pres- 
ence of magnetic fields [1-8]. The D3/D7 holographic 
system describes a confining 3-|-ld gauge theory -with 
quarks [9]. The magnetic field induces chiral symmetry 
breaking. The symmetry breaking and quark confine- 
ment are lost at high temperature and density. Bet-ween 
is a rich structure of phase transitions of both first and 
second order meeting at critical points. These transitions 
have been explored in [4] and the summary phase dia- 
gram is displayed in Fig la. Here the theory is interest- 
ing as a loose analogue for QCD -which is also a confining 
and chiral symmetry breaking gauge theory but -where 
-we can not as yet compute the precise phase diagram. 

Interest has also turned to the D3/D5 system [10] that 
describes fundamental representation matter fields on a 
2-|-ld defect -within a 3-t-ld gauge theory. This system 
may have some lessons for condensed matter systems. 
In [7] an analysis of the D3/D5 system at finite density 
(d) and at zero temperature (T) revealed that the chiral 
symmetry breaking transition -with increasing magnetic 
field {B) is not second order but similar to a Berezinskii- 
Kosterlitz-Thoulcss (BKT) transition [11] (see also the 
holographic example in [8, 12]). That is order parameters 
across the transition gro-w as exp(— a/^i/c — v) -where a 
is a constant and v = d/B. {vc the critical value for the 
transition). For small T the authors of [7] sho-wed the 
BKT transition returns to a second order nature. This 
difference from the D3/D7 case is surprising so it seems 
-worth fleshing out the entire phase diagram for the the- 
ory to see if other surprises are present. In this letter -we 
present that analysis - much of the computation matches 
that in the D3/D7 system -which -we -worked through in 
detail in [4] so here -we very briefly present the formalism 
and the conclusions. We display the resulting phase dia- 
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gram for massless matter fields in Fig lb. Clearly much 
of the structure is similar to the D3/D7 case but the 
second order boundary of the chiral symmetry breaking 
phase is distorted by the presence of the BKT transition. 



II. THE HOLOGRAPHIC DESCRIPTION 

The 3V=4 super Yang-Mills gauge theory at finite tem- 
perature has a holographic description in terms of an 
AdSs black hole geometry (-with N D3 branes at its 
core) [13]. The geometry can be -written as 



ds^ = -^{-gtdt^ + gxdx^ + gxdy^) 

+ ^(dp^ + p^dnl + dL^ + L'^d^l) , (1) 

-where x is t-wo dimensional, y -will be the D3 coordinate 
not shared by our D5, -we have split the transverse six 
plane into t-wo three planes each -with a radial coordinate 
p, L and a t-wo sphere, = AnggNa ^ and 



9t := 



2w* 



(2) 



The temperature of the theory is given by the position 
of the horizon, wh = ttR^T 

We include our 2-|-ld defect -with fundamental matter 
fields by placing a probe D5 brane in the D3 geometry. 
The probe limit corresponds to the quenched limit of the 
gauge theory. The D5 probe can be described by its DBI 
action 



S 



DBI — 



-Td5 J d'ev/-dct(P[G]a& + 2TTa'Fab 



(3) 



-where P[G]ab is the pullback of the metric and Fab is 
the gauge field living on the D5 -world volume. We will 
use Fab to introduce a constant magnetic field (eg F12 = 
—F21 = B) [1] and a chemical potential associated -with 
baryon number At{p) 7^ [14, 15] We embed the D5 
brane in the t, x, p and fl2 directions of the metric but to 
allo-w all possible embeddings must include a profile L{p) 
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FIG. 1: The phase diagrams for the D3/D7 [4] and D3/D5 systems, wh measure the temperature of the theory whilst fi is 
the chemical potential. The dashed line is a second order transition associated with the formation of quark density and meson 
melting. The dotted line is a second order transition for chiral symmetry restoration. In the D3/D5 case that transition ends 
at a BKT transition point and its effects on the second order line can be seen. The continuous line is the merged first order 
transition. The position of critical points are marked. 



at constant y,^2- The full DBI action we will consider 
is then 



dpL{p), (4) 



where £2 is a volume element on the 2-sphere and 



-NfT, 



W 



H 




(5) 



Since the action is independent of At , there is a conserved 
quantity d (^.~ ^^j^ and we can use the Legcndrc trans- 



formed action 
S^S- 
where 



dfFpt 



SS 

sfT 



pt 



{Is^^^ /t^^^a*) Jdp^ip) , (6) 



L : 
K 



-NfT, 



2^/2^^;4 



Kil + id,m (7) 



+ 



4 I 4 



d^ 



+ w%) {NfTob^TTa') 



l\2 



(8) 



To simplify the analysis we note that we can use the 
magnetic field value as the intrinsic scale of conformal 
symmetry breaking in the theory - that is we can rescale 



all quantities in (7) by B to give 



C = -NfTosiRVB) 
K^[^l±^]p^ 



'JL^K{l + U^), (9) 



1 



w 



~4 , ~4 P 



where the dimcnsionless variables are defined as 

iw,L,p,d) (11) 



L 



P 



RV2B' RV2B' {RVB)'^NfTD5'2TTa' J 

In all cases the embeddings become flat at large p tak- 
ing the form 



Hp) 



^ c 
m + — 

P 



(12) 



In the absence of temperature, magnetic field and density 
the regular embeddings are simply L(j)) = fh, which is 
the minimum length of a D3-D5 string, allowing us to 
identify it with the quark mass as shown, c should then 
be identified with the quark condensate. 

We will classify the D5 brane embeddings by their 
small p behavior. If the D5 brane touches the black 
hole horizon, we call it a black hole embedding, other- 
wise, we call it a Minkowski embedding. We have used 
Mathematica to solve the equations of motion for the D5 
embeddings resulting from (9). Typically in what fol- 
lows, we numerically shoot out from the black hole hori- 
zon (for black hole embeddings) or the p = axis (for 
Minkowski embeddings) with Neumann boundary condi- 
tion for a given d. Then by fitting the embedding function 
with (12) at large p we can read off fh and c. 

The Hamilton's equations from (6) arc dpd = and 



doAt 



5d ■ 



The first simply means that d is the con- 
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served quantity. The second reads as 



dpAt 



i + iuy 



-'H 



K 



(13) 



where A* 



V22-iTa'At 
R.V2B ' 



There is a trivial solution of (13) with d = and con- 
stant At [16]. The embeddings are then the same as 
those at zero chemical potential. For a finite d, A'^ is sin- 
gular at p = and requires a source. In other words the 
electric displacement must end on a charge source. The 
source is the end point of strings stretching between the 
D5 brane and the black hole horizon. The string tension 
pulls the D5 branes to the horizon resulting in black hole 
embeddings [14]. For such an embedding the chemical 
potential(pt) is defined as 




FIG. 2: A plot of the condensate vs the quark mass to show 
the first order phase transition at zero chemical potential in- 
duced by temperature. The solid fine corresponds to the black 
hole embedding and the dotted line to a Minkowski embed- 
ding. From bottom to top the curves correspond to temper- 
atures iuH = 0.25, 0.3435, 0.45. 



lim At{p) 

p~yoo 



dp d 



l\2 



PH 



w 



(14) 



H 



for a well defined At at the 



where we fixed At{pH) 
black hole horizon. 

The generic analysis below with massless quarks and 
B, T and /x all switched on involve four types of solution 
of the Euler Lagrange equations. All of these approach 
the p axis at large p to give a zero quark mass. Firstly, 
there are Minowski embeddings that avoid the black hole 
so have a non-zero condensate c - these solutions have 
d = so At = p. Secondly, there can be generic black 
hole solutions with both of c and d none zero. Finally 
there are solutions that lie entirely along the p axis so 
that c = but with d cither zero or non zero. In fact the 
flat embeddings with d ~ are always the energetically 
least preferred but the other three all play a part in the 
phase diagram of the theory. 

To compare these solutions we compute the relevant 
thermodynamic potentials. The Euclidcanized on shell 
bulk action can be interpreted as the thermodynamic po- 
tential of the boundary field theory. The Grand potential 
(51) is associated with the action (5) while the Helmholtz 
free energy (F) is associated with the Legendre trans- 
formed action (6): 

FiwH,d) 



NfTD5{RVB)^Yo\ 



oo — 4 — 4 
W — WJ 



dp 



^V(l + (^')^)Vi^ (15) 



Pff 



fl{wH,p.) := 



-S 



NfTDz{Ry/BfYo\ 



OO — 4 — 4 

- -7-;(i + (iO^)iii_::^(i6) 



where Vol denote the trivial 5-dimensional volume in- 
tegral except p space, so the thermodynamic potentials 
defined above are densities, strictly speaking. Since 
K ^ p'^, both integrals diverge as fP at infinity and need 
to be renormalized. 



III. CHIRAL SYMMETRY RESTORATION BY 
TEMPERATURE 



The chiral symmetry restoration transition by temper- 
ature is first order [6] (a transition related to the ther- 
mal transition for non-zero mass at B=0 [17]). The 
transition on the gravity side is between a Minkowski 
embedding that avoids the black hole to an embedding 
that lies along the p axis ending on the black hole. 
Fig 2 shows the (— c, m) diagram for some temperatures 
{wH = 0.25,0.3435,0.45 from the bottom). The solid 
lines are the black hole embeddings and the dotted lines 
are Minkowski embeddings. Since we are interested in 
the case to = 0, the condensate is the intersect of the 
curves with the vertical axis. As temperature goes up 
the condensate moves from the lower dot to the mid- 
dle curve continuously, then jumps at wh = 0.3435 to 
the origin (zero condensate), which corresponds to the 
chiral symmetric phase. It is also the transition from a 
Minkowski (dotted line) to a black hole embedding (solid 
line). This jump can be seen by a Maxwell construc- 
tion: TO and c are conjugate variables and the two areas 
between the middle curve and the axis are equal at the 
transition point. See [6] for more details. 

This transition as well as restoring chiral symmetry 
also corresponds to the melting of bound states of the 
defect quarks since the Minkowski embedding has sta- 
ble linearized mesonic fluctuation whilst the black hole 
embedding has a quasi- normal mode spectrum [18]. 
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(a) wh = 0.15 (b) wh = 10 ^ (c) wh = 10 ^ near the critical 

point 



FIG. 3: Plots of the condensate vs chemical potential on fixed temperature slices, showing the phase structure of the theory. 
Figure (b) and (c) show that at low temperature the BKT transition becomes second order. 



IV. CHIRAL SYMMETRY RESTORATION BY 
DENSITY 

At zero temperature we find two phase transitions with 
increasing chemical potential. 

At low chemical potentials the preferred embedding is 
a Minkowski embedding with = ^ so there is no quark 
density. There is then a transition to a black hole em- 
bedding with non-zero quark density, d. This transition, 
whilst appearing first order in terms of the brane embcd- 
dings, displays second order behaviour in all field theory 
quantities such as the condensate or density (which grows 
smoothly from zero). The transition also corresponds to 
the on set of bound state melting since the black hole 
embedding has quasi-normal modes rather than stable 
fluctuations. 

The chiral symmetry transition induced by density at 
zero temperature is distinct and also a continuous tran- 
sition. It has been shown to be of the BKT type for this 
D3/D5 case [7] as opposed to a mean- field type second 
order transition as seen in the D3/D7 case [4, 5]. 

The chiral symmetric phase corresponds to the trivial 
embedding, L = Q. Chiral symmetry breaking is signaled 
by the instability of small fluctuation around the L = Q 
embedding. The Free energy (15) with (9) at zero T 
reads 



F ~ ^l + L'^pi + ^+d? , (17) 
which can be expanded up to the quadratic order in L as 
F ^ A^l + pi + d^L'' + (18) 

At p ^ 1, 4 behaves as a scalar with = —2 in AdS4, 

while at small p <^ 1 and p <C d it behaves as a scalar 
with = ^ in AdS2 . The Breitenlohner-Freedman 

(BF) bound of AdS2 is -i, so below dc ^ V? the BF 
bound is violated and the embedding L = is unsta- 
ble [7]. This critical density corresponds to the critical 



chemical potential ~ 2.9 as can be computed from (14). 
In [7] it was shown that the condensate scales near this 
transition as 




which corresponds to BKT scaling [11]. This transition 
is an example of the analysis in [12] where it was shown 
that if a scalar mass in a holographic model could be 
tuned through the BF bound a BKT transition would be 
seen at the critical point. 

V. PHASE DIAGRAM IN fi-T PLANE 

To compute the full phase diagram we work on a series 
of constant T slices. We have found the four relevant em- 
beddings discussed above and found those that minimize 
the relevant thermodynamic potential. For more details 
of the method and relevant analysis we refer to [4] , where 
we studied D3/D7 system using the same methods. Fig 3 
shows some example plots of the dependence of the con- 
densate on the density on fixed T slices. It shows that the 
Minkowski embedding with d ~ is preferred at low /i, 
a black hole embedding with growing d at intermediate 
/2, before finally a transition to a flat embedding occurs 
at high chemical potential. 

Qualitatively the phase diagram, shown in Fig 1, is 
almost the same as the D3/D7 case - the two second 
order transitions at zero temperature converge at two 
critical points to form the first order transition identi- 
fied at zero density. The only difference is induced by 
the chiral phase transition at zero T. Comparing to the 
D3/D7 case we see there is a long tail near zero T, the end 
point of which corresponds to the BKT transition. How- 
ever even infinitesimal temperature turns it into mean- 
field type second order transition[7, 8]. In Fig 3bc we 
plot the condensate against p at a very low temperature 
{wh = 10~^)to show the second order nature. 
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